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The independent Hall, magnetoresistive, and thermoelectric coefficients for a tetragonal 
crystal have been tabulated and geometric configurations for their experimental measure- 
ment have been determined. These coefficients have been calculated on assumptions of 
several simple ellipsoidal models, in the range of nondegenerate statistics. Implications 
of experimentally observed isotropy or anisotropy of transport properties on the structure 
of the energy surfaces are noted. 



1. Introduction 

Detailed studies on the transport properties of 
rutile require tabulations of independent tensor 
coefficients. In this report we display these trans- 
port coefficients and method of measurement for 
galvanomagnetic and thermoelectric effects in tetrag- 
onal crystal structures. In anticipation of future 
experimental and theoretical results these transport 
coefficients have been calculated on the basis of 
various ellipsoidal energy band structures compatible 
with tetragonal symmetry. 



2. Galvanomagnetic Coefficients 

The electric current density components J t and 
electric field components Ei are related by the 
conductivity tensor a tj in the following way: 



J,=cr„(H)E„ 



(1) 



where the conductivity is understood to be a function 
of the magnetic field vector H and the Einstein sum- 
mation convention is used. For weak field galvano- 
magnetic effects a series expansion of a in powers of 
H up to second order is frequently sufficient. To 
second order the current is given by 



J t =a<$E J +*l%E j H k +*l% l E J H k H l . 



(2) 



The superscripts indicate the corresponding power of 
H, but these will frequently be omitted without 
causing ambiguity. The restrictions of crystal 
symmetry [l] 1 and the Onsager reciprocal relations 
[2, 3] greatly limit the number of independent con- 
ductivity coefficients in eq (2). The symmetry 
conditions of the tetragonal point group of rutile 
(D Ah ) require that eq (2) remain invariant under 
rotations about three perpendicular 2-fold axes, one 



[ Figures in brackets indicate the literature references at the end of this paper. 



of which is also a 4-fold axis. The Onsager rela- 
tions require that alf and <r$ kl be symmetric, while 
<r\)\ must be antisymmetric on interchange of i and j. 
The resulting independent components are produced 
in tables 1, 2 and 3. Note that for zero order there 
are two coefficients, for first order two, and for 
second order there are seven. 



Table 1. Zero order magneto-conductive coefficients 

Asterisks denote the coefficients which will be taken as the independent quanti- 
ties. Indices 1, 2, 3 refer to crystal axes, 3 being the 4-fold axis. 



All others vanish. 



Table 2. First order magneto-conductive coefficients 



All others vanish. 
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Table 3. Second order magneto-conductive coefficients 
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cr 3311 ==0 '3322 



All others vanish. 



More often than not, experimental measurements 
yield resistivity coefficients, the reciprocal of the 
conductivity. The resistivity coefficients are given 
up to the second order by 



E i =[ P if+ P l)lH k +p%H k H t ]J J . 



(3) 
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The same symmetry conditions apply as before, and 
the independent resistive coefficients will be identified 
by the same subscripts. The coefficients are given 
in table 4. 

Table 4. Magnetoresistive coefficients in terms of magneto- 
conductive coefficients 

Expressions are correct to second order in the magnetic field 
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3. Experimental Determination of Magneto- 
resistive Coefficients 

The zero order resistivities are measured in the 
usual way. The two coefficients p n and p 33 refer to 
the resistivities perpendicular and parallel to the 
4-fold axis, respectively. Typical experimental ori- 



entations and geometries have been determined for 
obtaining the magnetoresistive coefficients. The 
first order, or Hall, coefficients are described in table 
5, and the second order in table 6. Numerals 1, 2, 3 
in the tables refer to crystal axes; 1 and 2 being 2-fold 
and 3 being the 4-fold axis. 



Table 5. Configurations for measurement of Hall effect in 
tetragonal cry 
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Table 6. Configurations for measurement of magnetoresistive 
coefficients of tetragonal crystals 
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4. Magnetoconductive Coefficients on the 
Basis of Simple Ellipsoidal Models 

Experimental and theoretical studies of band 
si ructure and transport coefficients may be aided by 
calculations of the transport properties of simple 
many ellipsoid energy surface configurations. This 
type of analysis of magnetoresistance has been very 
successful in the determination of the band structures 
of silicon and germanium [4, 5]. 

First we discuss the magnetoconductive coefficients 
of a single ellipsoid. Rather than using a Boltzmann 
equation, we follow Brooks [5] method of following 
the motion of a single charge carrier on the energy 
surface under the influence of applied fields, and 
then averaging over a Maxwell distribution of 
carriers. The equations of motion for an electron 
on an ellipsoidal energy surface, in the effective 
mass approximation, are the three equations 

iK^+Wc^-^' (4) 

where the subscripts indicate components or effective 
masses along principal directions of the ellipsoid, 
and i, j, k are cyclic permutations of 1, 2, 3. These 
equations of motion are solved to second order in 
H, giving the velocity as a function of time. The 
acquired motion is obtained by averaging over a 
distribution of free drift times given by (\Jr)e~ t,T . 
The resultant acquired current averaged over the 
Maxwell distribution of electrons is 



m^ miTrijC mi7n k c 



V) 



mim H c 



E t H* r 



m\m,jc 2 



E f Hl 



mitrijm k c 



2 E k H t H k + e ' {rZ) v EjHjH*. (5) 



In eq (5) it is understood that the relaxation time r 
is independent of direction, but perhaps dependent 
on energy. The angular brackets indicate averages 
of r over the Maxwellian energy distribution. The 
total current is just the average current multiplied 
by the number of electrons in the ellipsoid. 

We shall consider now the resistivity coefficients 
for several configurations of ellipsoids compatible 
with crystal symmetry. The assumption is made 
that the total conductivity of the array is the sum 
of the conductivities of the individual ellipsoids. 
Thus we assume that a relaxation time exists for 
each ellipsoid, and that intervalley scattering is 
negligible. 



4.1. Single Ellipsoid 

In this case mi^^7n 2 9 £ m^ The configuration is 
shown in figure 1. The results for the resistive 
coefficients are as follows: 




m l z m 2 ^m 3 



Figure 1. Single ellipsoid for tetragonal symmetry. 
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4.2. Four Ellipsoid Model 

In this case the ellipsoids can have three different 
masses, but are placed along the 2-fold axes, as 
shown in figure 2. The effective masses m 1; ra 2 , ra 3 
apply to ellipsoid number 1 of figure 2. The masses 
for the other ellipsoids are obtained by rotation. 
The resistive coefficients are: 
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Figure 2. Four ellipsoids for tetragonal symmetry. 
4.3. Five Ellipsoids 

There is a possibility that more than one energy 
surface will be populated. To give an example, we 
shall consider a combination of cases (a) and (b), 
as shown in figure 3. The resultant resistivities are: 
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Figure 3. Five ellipsoids for tetragonal symmetry. 
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In all the previous formulas the temperature 
dependence of the resistivities depends upon averages 
of several powers of r. For power law scattering, 
T=r t~ p , where e is the energy and p the power, we 
have 



(r»>= 



(kT)-'">T (^-np) 

r (D 



(6) 



It might be worth noting that in tetragonal crystals 
anisotropic effects make themselves known in lower 
order magnetoresistive coefficients than for cubic 
cases. Thus, for example, information concerning 
anisotropy of the energy surfaces can be obtained 
from the Hall coefficient of Ti0 2 , while one must go 
to magnetoresistive to see the lowest order effects 
of cubic anisotropy in Si or Ge. 

5. Thermoelectric Power 

The electric current density components J t and 
the thermal current components Q t are related to the 
electric field components E t and the temperature 

dT . 
gradient components r-^r in the following way: 



Ji—AYjEj-^A™ 
Q i =Af j E j J r Af j 



dT 
bT 

dx; 



(7) 



(8) 



However, in order to apply Onsager's theorem to the 
above coefficients we mnst rewrite the equations in 
terms of affinities and fluxes [6]. Thus (7) and (8) 
become 



_Jj = R11 I i>M R12 _d_ 1 



Qt=m i , 
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d 1 



11 T dXj ' " bXj T 



(9) 
(10) 



where e is the charge of an electron and n is the elec- 
trochemical potential. The Onsager theorem states 
that B\)=B l *i and thereby reduces the number of 
independent coefficients. Accordingly, the number 
of the independent coefficients A\) is also reduced. 
Applying the restrictions of crystal symmetry to 
the coefficients A\) it follows that there are two 
independent elements A u l =AU and A\\. 

The absolute thermoelectric power is defined as 
the change in voltage per unit change in temperature 
difference at zero current. So setting J t equal to 
zero and using the above results we get: 



E r - 



A\* t dT 

"a\\ dx; 



an 



The voltage V is given by 



Vi =j ti E l( I Xr - -^^^.^ 



'o A 12 

u IT 
An ai i* 



The absolute thermoelectric power is 

A\T A12 

u ' v i -^ii or 

dT~ A\\~ " 



(12) 



(13) 



It follows that we have two possible values of the 

A 1 ' 



A\t 



corre- 



thermoelectric power, ~-in an d 

sponding to the one and three axes respectively. 
We see also that the Onsager theorem does not give 
us any limitation on the possible numbers of values 
of the thermoelectric power but only gives us rela- 
tions between it and other thermoelectric coefficients. 



6. Thermoelectric Power on the Basis of 
Simple Ellipsoidal Models 



[71: 



We begin with the Boltzmann transport equation 



% ^ •V*J+? ) r V r(e) 



(14) 



Since the electric field and temperature gradient are 
always small and their squares and products can be 
neglected, we obtain an approximate solution of 
eq (14) by putting /==/ on the left side where f is 
the equilibrium distribution. Using 



»/o_ r l /e-r\ 5/q 
)T bT\ T ) de 



Wo 
dT 



we get 



f=f>+-" j si(mi+ T mw(T) + ' E >) (15 



5) 



where f is the Fermi energy, and e the energy. 
The electric current density is given by [7]: 



Jz= 4^ yJ dh 



(16) 



Using our solution for / and writing in terms of 
components of J we get: 

^= k ^ 2e ^ t Mt)1 +l SS} (i7) 
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where 



rv au -^— dk- 



47T 3 Jo 
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(18) 



(19) 



Setting </*=(), we calculate the thermoelectric power 
in the same manner as section 5. The result is: 



f ~ el K u }T 



(20) 



(Here we anticipated the results K ij =L i:j =0 if i^j.) 
Assuming a single ellipsoidal energy surface and 

using dk=- — A — as the element of volume where 
grad^e 

dS is an element of surface of the ellipsoid, w^e find 
upon performing the indicated integration: 
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(21) 



We [have again assumed t=t €~ p . (The above 
results are obvious on comparison with the expression 
for a a by Abeles and Meiboom [3].) The thermo- 
electric power for a single ellipsoid ma} r be written: 
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(22; 



The case for many ellipsoids is easily gotten from 
the above result. Since conductivities simply add 
we have the result: 

k 

where k refers to the kik ellipsoid. We note that 
if f, measured from the bottom of the ellipsoid, 
is constant for all ellipsoids of the configuration 
then the thermoelectric power is independent of 
direction and is given by 



S-- 



-h[(l-p) kT -^} 



(24) 



Thus the simplest case which could yield an aniso- 
tropic thermoelectric power and also have tetragonal 
symmetry would be the five ellipsoid case discussed 
earlier. 

Using the relationship 



6 r/*7=l 7^3(2^^-3/2 



(25) 



we can express the thermoelectric power in terms 
of the density of electrons, temperature, and effec- 
tive masses. 

6.1. Single-Ellipsoid Case 

The masses are those indicated in figure 1. Using 
eqs (24) and (25) we get for the case of a single 
ellipsoid: 



S i =S : ^ — G? 



where 



G=2k-k log Ti+I k log T-k log ^K 6 {2irmk)-^ 2 



and 



m= (m 1 m 2 mo) 1/3 . 
6.2. Four Ellipsoid Case, Figure 2 



(26) 

(27) 
(28) 



Carrying out a similar calculation for the four 
ellipsoid case we get: 

G=2k-k log |+| k log T-k log ~ h\27rrnk)- 3/2 , 

with eqs (26) and (28) still true. Note that the 
only change is the division of the density of electrons 
by 4. This is just what we would expect since 
equal voltage sources in parallel produce a voltage 
equal to one of the sources. Thus the presence of 
many ellipsoids is seen as an effective reduction 
of the number of carriers. 

6.3. Five Ellipsoid Models 

We shall consider the simplest anisotropic case of 
five ellipsoids as shown in figure 3. Using eqs (23) 
and (25), we obtain for the two values of the thermo- 
electric power: 
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1 



G l =2k-k log ^+| k log T-k log i h 3 (2Tm 1 k)- 2 

G 2 =2k-k log n,+§ k log T-k log | K\2imi 2 k)-v 2 

m 1 =(m 1 m 2 m 3 ) 1/3 
m 2 =(mlm 5 ) 1/3 . 
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c 

Here S a = is the thermoelectric power perpen- 

C 
dicular to the 4-fold axis and is the thermoelec- 

e 

trie power along the 4-fold axis. 

7. Phonon-Drag Contribution to the Thermo- 
electric Power 

In all the former work, the phonon-drag contribu- 
tion [8, 9] to the thermoelectric power was ignored. 
According to Frederikse [9] the thermoelectric power 
is given by: 



wli 



Iff 

'e'Ki 



is the extra term due to plionon drag. 



For the case of lattice scattering (2> =1 A), i.e., 



K_ 
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1 3 k vh (2ms>y 
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s= speed of sound, 
A P a=phonon-ph.onon mean free path, 
A cv =electron-phonon mean free path, 
A ei = electron-impurity ion mean free path, 

m=(m 1 m 2 m 3 ) 1/3 . 

We have the result that the thermoelectric power of 
a single ellipsoid is given by: 



Any easy extension to the case of many ellipsoids is 
obtained by use of [10] 



2W 

o _ 

total -4-f l 



(32) 



Here S k indicates the thermoelectric power of 
a single ellipsoid and <r\ its conductivity in the ith 
direction. This is easily seen to be equivalent to 
our former method of treating the case of many 
ellipsoids. 



Note: The formulas of sections 4.3 and 6.3 are 
derived on the basis of a single relaxation time for 
the different bands. Generalization to different 
relaxation times for each band is straightforward. 
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